Qubit protection in nuclear-spin quantum dot memories 
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We present a mechanism to protect quantum information stored in an ensemble of nuclear spins 
in a semiconductor quantum dot. When the dot is charged the nuclei interact with the spin of the 
excess electron through the hyperfine coupling. If this coupling is made off-resonant it leads to an 
energy gap between the collective storage states and all other states. We show that the energy gap 
protects the quantum memory from local spin-flip and spin-dephasing noise. Effects of non-perfect 
initial spin polarization and inhomogeneous hyperfine coupling are discussed. 



An essential ingredient for quantum computation and 
long-distance quantum communication is a reliable quan- 
tum memory. Nuclear spins in semiconductor nanostruc- 
tures are excellent candidates for this task. With a mag- 
neton 3 orders of magnitude weaker than electron spins, 
they are largely decoupled from their environment, and 
the hyperfine interaction with electron spins allows one 
to access ensembles of nuclear spins in a controlled way 
P, H H II H @, 0, B Hi • In particular, the quantum state 
of an electron spin can be mapped onto the nuclear spins, 
giving rise to a long-term memory 0, 0, 0, 0] ■ Neverthe- 
less, memory lifetimes are limited, e.g., by dipole-dipole 
interactions among the nuclei. In this Letter we demon- 
strate that the presence of the electron spin in the quan- 
tum dot substantially reduces the decoherence of this 
collective memory associated with surrounding nuclear 
spins. The virtual transitions between electronic and 
nuclear states can be used to produce an energy shift 
proportional to the number of excitations in the storage 
spin-wave mode. This isolates the storage states energet- 
ically and protects them against nuclear spin flips and 
spin diffusion. 

Consider a quantum dot charged with a single excess 
electron as indicated in Fig. [TJ The electron spin S is 
coupled to the ensemble of underlying nuclear spins V 
by the Fermi contact interaction, 
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where A is the average hyperfine interaction constant, 
A ~ 90 /ieV for GaAs, and Qj is proportional to the elec- 
tron density at the position of the jth nucleus, ^ . Qj = 1. 
For convenience, we introduce the collective operators 
A = Y^,j £?jl j • The first term in Eq. ([TJ) provides an effec- 
tive magnetic field B® H = A(A Z ) j ' g* hb for the electron, 
known as the Overhauser field. The same also produces 
an energy shift for each nuclei, the so-called Knight shift. 
The flip-flop terms in Eq. Q), H ic = ^(A+S-+A-S+), 
can be used to polarize the nuclear spins [J, 0] , and to 
map the electron's spin state into a collective spin mode 
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FIG. I: (Color online) Left: Charged quantum dot with a 
single, polarized excess electron. Right: Spectrum of the ef- 
fective nuclear Hamiltonian in the presence of a polarized elec- 
tron. Off-resonant hyperfine coupling results in a gap A gap 
between the storage state |1) and the non-storage states 
Ak denotes the Zeeman shift due to the effective magnetic 
field associated with the electron spin (Knight shift). 



of the nuclei [3|,|4j. As will be shown here, the same can 
be used to provide a protective energy gap. 

Fully polarized nuclei. We start by reconsidering the 
storage of a qubit in a collective nuclear state In 
the simplest case when all the nuclear spins are initially 
polarized in the — z direction (zero temperature limit), 
the |J.) e and ||) e spin states of the electron are mapped 
onto the nuclear spin states 
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respectively. iJjc couples the state |0)|f) e to |l)||) e with 
an angular frequency fl — A{J2j g^/) 1 / 2 . The detuning 
between these two states, S = S cl + (5 OH , comes from the 
electron's intrinsic energy splitting S el due to, e.g., an 
external magnetic field, and from the Overhauser field, 
S oli = —AI. Coherent flip-flops between the electron 
and nuclear spins can be brought into resonance [5 -C 0) 
through d) el , e.g., applying a spin-state dependent Stark 
laser pulse [10j]. Then |0)(a||) e +/3|T) e ) 

can be rotated 

to (a|0)+/3|l))|J,) e , and the quantum information can be 
transferred from the electron to the nuclear spin ensemble 
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and back 0,||. 

Assume that, after the qubit has been written into the 
nuclei, the polarized electron is not removed from the 
dot but the hyperfme flip-flops are tuned off-resonant 
(8 ^> Q). Now real transitions can no longer take place 
between |l)||) e and |0)|1%. However, the residual vir- 
tual transitions repel the two states from each other, in 
analogy to the dynamic Stark effect. As a result, af- 
ter eliminating the electron, the energy of state |1) gets 
shifted by A gap = — f2 2 /4i5. The other, orthogonal states 
also having exactly one spin flipped (denoted by \l q ) in 
Fig. [T]) are "subradiant" , i.e., are not coupled via Hjc to 
the electron. Therefore, they are unaffected by the shift. 
This is the origin of the energy gap. 

To understand the protection scheme, let us introduce 
nuclear spin waves. As long as the nuclei remain highly 
polarized, one can introduce bosonic operators through 
the Holstein-Primakoff transformation: hj w P__j\f^J, 
a}- « P + /y/2I, and a^aj = P z + I. This allows us to 
define the bosonic spin waves 
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where the unitary matrix r\ q] - describes the mode func- 
tions. We identify the storage mode q — as the one 
given by 770^ = V2I^gj, and write |1) = $J|0). This is 
the mode which is directly coupled to the electron spin. 
In fact, Hjc ~ §(^0^- + ^o<S+) 1S a Jaynes-Cummings 
coupling in the bosonic approximation. After eliminat- 
ing the electron, Hjc reduces to H ga , P — — =jgA + A- « 
Agapl'o'&o- As shown in Fig.[]] H gap lifts the degeneracy 
between states of different number of storage-mode exci- 
tations. This is the key feature of our protection scheme: 
any decoherence process that is associated with a tran- 
sition from the storage mode $0 to any other mode $ 9 
now has to bridge an energy difference. If this gap is 
larger than the spectral width of the noise, the effect of 
the noise is substantially reduced. 

A more detailed analysis shows that the off-resonant 
interaction with the electron spin — which itself is cou- 
pled, e.g., to phonons — leads in general also to an addi- 
tional decoherence mechanism for the nuclear spins. If 
the corresponding electron spin dephasing rate 7 is small 
compared to the electron's precession frequency 8, the de- 
cay rate for the storage mode is reduced by the low prob- 
ability of exciting the electron spin state: 7f2 2 /S 2 <C 7. 

In addition to the gap, the electron is also respon- 
sible for the Knight shift H K = AA Z (S Z ). The dif- 
ference of the Knight shifts for the |0) and |1) states, 
A K = -^Y.jffj/Y.jQ^ is typically much less than 
Agap- When the hyperfme coupling is inhomogeneous, 
however, |1) fails to be eigenstate of the Knight shift 
Hamiltonian: H K \1) = (-±<S OH + A K )|1> + C|1 X )> where 
the state {l- 1 } is orthonormal to |1) and the coupling pa- 
rameter C 2 = 4 J2j 6j I J2j o] ~ Aj^ characterizes the 



inhomogeneities. As a consequence, the storage mode is 
only an approximate eigenmode, and it gradually mixes 
with non-storage modes as time passes. This causes loss 
of the stored qubit. If -1 ) is, however, off-resonant due 
to the energy gap, and our simulations show that the 
corresponding probability of finding the system in state 
|l x ) is bounded by 4C 2 /A 2 ap , so the detrimental effect 
of the inhomogeneous Knight shift is suppressed by the 
energy gap. In addition, since the admixture of jf" 1 ") is a 
coherent process, it can be cancelled by refocusing (echo) 
methods. 

A large gap can be achieved by bringing the hyperfme 
interaction close to resonance. For example, a non-zero 
external magnetic field or laser induced AC Stark shifts 
flpj ] can partially cancel the Overhauser field, such that 
8 <C 5 C \ ps — <5oh = AI. (Of course, 8 should be kept suf- 
ficiently large so that the hyperfme coupling remains off- 
resonant). The requirement of separation of time scales 
implies C < | A gap | -C fi -C \S\, i.e., 8 > 10Q. To estimate 
the orders of magnitude of the different energies, we take 
an oblate Gaussian electron density of ratio (1,1,1/3), 
and we consider spin-^ nuclei. Then it is easy to see that 
Ak and Q are inversely proportional to the number of 
nuclei TV, whereas VI, A gap cx TV -1 / 2 only (Fig. [2^). 

To analyze the decoherence suppression, we first con- 
sider a simplistic noise model where the nuclear spins are 
coupled to fluctuating, classical fields. The correspond- 
ing interaction Hamiltonian is given by V — Y] j B J • P . 
We assume isotropic Gaussian noise with zero mean and 

Bl(t)BW) = 5 MV tikCe-r^ (5) 

for fi, v = x,y, z, where specifies the spatial correla- 
tions of the noise acting on different nuclei. For simplic- 

Linear size of dot (a) [in lattice unity] 
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FIG. 2: Hyperfme Rabi frequency (fi), protective energy 
gap (Ag a p), Knight shift difference between the logical states 
(Ak), symmetry breaking couplings due to inhomogeneities 
(£ and uj), qubit decoherence rate due to dipolar spin diffu- 
sion without (Td) and with (r^) protection, (a) The fully 
polarized (zero temperature) case is displayed as function of 
the number of spin-i nuclei (TV) taking part in the storage, 
i.e., located within 3ct of the oblate Gaussian electron distri- 
bution with in-plane variance a. (b) Estimated energies in 
dark states |©n,/3) with n spins flipped from the fully polar- 
ized state for ./V = 10°. Energy units are chosen to match 
GaAs. 
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N = 10 3 > 



TV = 10 4 / 
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Noise correlation lenght (£) [lattice units] 

FIG. 3: The parameters H and S' describing the effects of 
spatial correlations in the classical noise = e~ Tjk ^) for 
different number of nuclei. The same family of Gaussian elec- 
tron densities was used as in Fig. The bullets on the curves 
denote the linear size of the dot given by the variance a. 



ity, the noise spectrum is assumed to be Lorentzian with 
a width r, although similar results hold for other spectra 
with a high-frequency cut-off. 

Let us first discuss the dephasing part, V z — B Z P Z , 
of the noise. Using the bosonic spin-wave operators in- 
troduced in Eq. Q we can express V z as 



V z = 



Vpj Vqj 



(6) 



Dephasing of individual nuclear spins thus means transfer 
of excitations between different spin-wave modes. Espe- 
cially, it leads to both real and virtual transitions from 
|1) to a non-storage state (with q ^ 0). As the latter 
state is "subradiant" and, thus, equivalent to |0) when 
the memory is read out, this process essentially results in 
damping (for real transitions) and dephasing (for virtual 
transitions) of the stored logical qubit This can be 
seen by formally eliminating the classical fields and all 
non-storage mode in Markov approximation and deriv- 
ing a master equation for the storage mode. For that, we 
assume the zero temperature limit with all non-storage 
modes $ 9 /o m the vacuum state. This results in 



dt 



P = i 



(7) 



with energy shift E z = (1 - S)CA gap /(r 2 + A 2 ap ) and 

£z [p] = 7i (2$oP$o ~ $o*oP - /5$o*o) 
+ 72 (2*J$ Mo$o ~ H$oH$oP ~ p^oH^o) ■ (8) 

Here, 71 is the damping rate of the stored qubit while 72 
describes its dephasing. The two rates are given by 



CT 



7i 



r 2 + a 2 

' gap 



(1-3), 
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where we have introduced the dimensionless parameter 
s = Y, 3 k^okQ 2 jQ 2 k l '{YjI BiY containing the spatial part 
of the noise correlator. 



When the correlation length of the classical noise is 
smaller than the distance between the nuclei (local un- 
corrected noise, ~ Sjk), S scales inversely with the 
number of nuclei (Fig. [3]). In this case, the dephasing 
rate 72 vanishes as 1/N, which is an effect of the collec- 



tive nature of the storage states The storage of a 

qubit corresponds to an encoding of the logical state in 
a large, delocalized ensemble of N physical spins. As the 
decoherence has strongly local character, there is only a 
very small effect on the dephasing of the qubit. Secondly, 
the loss of the stored qubit is due to transitions among 
states with different number of excitations in the storage 
mode. These transitions are strongly suppressed and the 
damping rate 71 is decreased if A gap is large compared 
to the width of the noise spectrum T (or the correspond- 
ing cut-off frequency) . Finally, we note that the opposite 
limit of infinite spatial correlation length (^jk = 1) cor- 
responds to a homogeneous random field resulting, e.g., 
from a global external source. In that case, as 1 (see 
Fig. [3]) and there is no protection against dephasing. 

Following a similar but slightly more involved pro- 
cedure we can discuss the spin-flip part V xy = 
{B J + I J _ + B 3 _I J + ^ of the noise. When deriving a 
master equation for this case, we need to keep higher 
order terms in the Holstein-Primakoff approximation: in 
the next order P_ w v2j(l — Xc^-aj)a,j (and similarly 
for P + ) with A = 1 - (1 - 1/2/) 1 / 2 . Here we have ne- 
glected the probability of double or more excitations on 
the same site j, which is reasonable in the high polar- 
ization (T = 0) limit and exact for spin-1 nuclei. Omit- 
ting the energy shifts, the Lindbladian describing deco- 
herences due to spin flips reads, in leading order of 1/N, 

£xy(p) = (73 + 74)(2<fro/5&o - H$oP - pH$o) 
+ 7 5 (2$J$ p$ ( ! ) $o - $i$oH$oP ~ p&o&oH&o) 

+ 7 3 (2$Jp$o-MSp-pMS), (10) 

which describes decay with rate 74, dephasing with rate 
75, and additionally thcrmalization (relaxation to the 
identity matrix) with rate 73. The rates read 
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2CT/A 2 



(A g 



75 



r 2 + A|(E^: 



A K ) 2 ' 
(11) 



In the limit of vanishing spatial correlations of the spin- 
flip noise, 5' = J2jk^kQjQk/Y,i Qi tends to 1 (Fig-© 
and we have protection against thermalization (73) be- 



cause of the separation of |0) and |1) by A { 



ap 



The 



decay corresponding to 74 is due to spin-flip induced tran- 
sitions between |1) and |l p , l g ) (the latter containing a 
total of two excitations but none in the storage mode), 
and the energy to bridge is in the order of A gap — Ak 
(see Fig.[T]). Finally, the last factor in the dephasing rate 
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75 scales as 1/N, indicating that it is the collective na- 
ture of the storage that leads to protection. Note that 
the nonlinearity of the Holstcin-Primakoff representation 
is responsible for this dephasing: the virtual non-storage 
excitations are interacting with the storage mode. 

Another potential source of decoherence is nuclear spin 
diffusion due to dipole-dipole interaction between nuclear 
spins [12j. The energy gap gives protection against this 
effect, too. The dipolar interaction between the pairs of 
spins is described in the secular approximation by 



H 
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B 



A 1 - 



2I^B jk a]a k , (12) 



where Bj k = j"f 2 {3 cos 2 9j k — l)/?*| fe , 7 is the gyromag- 
netic factor, = rj — r k is the distance between two nu- 
clei, Ojk is the zenith angle of the vector Tj k , and we used 
the first order Holstein-Primakoff approximation. At full 
polarization, we can rewrite the dipolar Hamiltonian Q12p 
in terms of the bosonic spin wave mode operators ([4]) 

as H D = J2pq B pq^q, with B P q = J2j=ik B jkVpjVq k - 

Thus, the storage mode is coupled to a bath of non- 
storage modes as if it were a central spin coupled to 
a mesoscopic spin bath [l3[ 14 1. Although the total 



number of excitations is conserved, Hjj is responsible 
for decoherence of the qubit via transitions from the 
storage state |1) to non-storage states In fact, 

the non-storage modes produce a fluctuating effective 
transversal magnetic field with (complex) Larmor fre- 
quency f2ff_ = 2 J2 q ^o B 0q^q- If the electron were not 
present, these fluctuations would lead to a decoherence 
rate T D ~ AO^ ff = {2J2 q ^o \B 0q \ 2 ) 1/2 in the fully po- 
larized state, which is numerically found to be in the 
order of 100 Hz for GaAs (Fig. gji). With the protec- 
tive gap, however, the storage mode creation and anni- 
hilation operators ($J and $0) rotate rapidly with re- 
spect to the other ones, and the above coupling aver- 
ages out and disappears in first order of the dipolar 
perturbation. In second order, the strength of the re- 
maining coupling between the storage mode and mode q 
is proportional to Aj ap X^r^o B 0rB rq , and the corre- 
sponding fluctuations yield a decoherence rate of T' D ~ 
A ga 1 p (2E^o|E, # o^S rg | 2 ) 1 /2 „ 3 x 10W/A gap . 
Typically, A gap ~ 1 MHz depending on the dot size 
(Fig. [2^), so the effects of spin diffusion can be suppressed 
by several orders of magnitude. 

Non-perfect spin polarization. Finally, we investigate 
the consequences of non-perfect nuclear spin polariza- 
tion. It has been shown that partially polarized nu- 
clei (at finite temperature) can also be used for storing 
a qubit state [J]. Instead of the fully polarized state 
@, the initial preparation drives the nuclear ensemble 
into a statistical mixture of dark states \D n-0 ) defined by 
A-\T> nt p) = 0. These dark states can be characterized 
by the total number of spins flipped n and the permu- 
tation group quantum number f3. As the detuning 5 is 



adiabatically swept from far negative to far positive, a 
superposition of the ||) e and ||) e electron spin states is 
mapped into the mixture of superpositions of the nuclear 
spin states \D n ,0) an d \£n,p) = jF^+l^n./j), and the 
qubit state is efficiently written into the memory [JJ. 

When the electron is left in the quantum dot, it feels 
different Overhauser fields for different dark states, hence 
the detuning should be adjusted such that S° H + S cl » 
Var(<5„ ). Moreover, the hyperfine Rabi frequency also 



varies with n and the energy gap A g 



is not the 



same for all dark states. This inhomogeneous broaden- 
ing would result in dephasing of the qubit, but can be 
avoided by a symmetric spin echo sequence Q. 

To describe inhomogeneous effects in the case of non- 
perfect polarization, first we note that the storage state 
\D n ,p) is no longer an eigenstate of the Knight shift oper- 
ator, but it is partially mapped into an orthogonal state: 
H K \T>„,p) = -\8° H \V nt0 )+u n \V^ ). This is due to the 

fact that the inhomogeneous A Zj ± operators do not follow 
the angular momentum commutation relation. Further- 
more, \£ n ,p) is neither an eigenstate of -ff gap nor of Hk'- 
H\£ n ,p) = (-±<5° H + A K , n + A gap ,„)|£„ i/3 ) + ( n \£^ ). 
The parameters can be expressed as expectation values 
in \T> n: p): 

n 2 n =A 2 (A_A + ), u; 2 n = ^((A 2 z )-(A z ) 2 ), 

= A i {A-A+A-A+) /46 n n 2 n , 

A K ,„ = 4(A 2 ) - A 3 {A_A z A + }/2Q 2 n , 

( 2 = {£ n ,p\H 2 \£ n ,p) - (£ nt0 \H\£ n>0 ) 2 . (13) 

The explicit form of the inhomogeneous dark states [3| 
allows us to estimate these values (see Fig. [5b). We ex- 
pect that the storage mode is still protected as long as 
to n and are much smaller than A gap , which is the case 
even for considerable unpolarized fraction (n/N). 

In summary, we have demonstrated that it is possible 
to suppress the influence of spin-dephasing and spin-flips 
on a quantum memory consisting of a delocalized ensem- 
ble of nuclear spins in a quantum dot if the noise has 
a highly local character and the spectral width or cut- 
off frequency of the noise spectrum is small compared to 
the energy gap. We have shown in particular that the 
memory can be protected against nuclear spin diffusion 
mediated by dipole-dipole interaction. We have also an- 
alyzed the effects of inhomogeneous hyperfine couplings 
and imperfect initial nuclear spin polarization. 
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